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We investigate the infrared behaviour of the gluon propagator in Landau gauge on a lattice with 
twisted boundary conditions. 

Analytic calculations using Dyson-Schwinger equations, exact renormalization group and 
stochastic quantization show that the gluon propagator in Landau gauge approaches zero for small 
momentum. On the other hand lattice calculations and calculations on a four-torus seem to give 
rise to a non-zero limit. One possible reason for this difference is the existence of zero-momentum 
fluctuation modes which potentially give a massive contribution to the gluon propagator. Our sim- 
ulations show that with twisted boundary conditions these zero-momentum modes are suppressed 
and the gluon propagator becomes smaller than in a periodic ensemble. 
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1. Introduction 

On the gluon propagator has been calculated with several techniques, using Dyson-Schwinger 
equations [[I]], stochastic quantization and exact renormalization group equation [||]. Represent- 
ing the gluon propagator in Landau gauge as 

Df,{p) = 5"\5^,-p,p,/p^)Dip^) (1.1) 

the low momentum behaviour was found to be D{p^) ''oc" (p2)(2K--i) ^j^j^ ^ ^ other 
hand calculations on the torus using Dyson-Schwinger equations [Q] and lattice calculations 
[Q, ^, ^, ^, ^ seem to indicate a non-zero value of D at p = 0, i.e. K = 0.5. A possible explanation 
for the difference between R'^ (non-compact space-time) and (compact space-time) are zero 
momentum modes on the compact space-time manifold which potentially give rise to a massive 
contribution to the gluon propagator. 

In this talk we wish to investigate how the implementation of non-trivial boundary conditions 
changes the infrared behaviour of the gluon propagator. 

2. Twisted boundary conditions 

To determine the influence of different boundary conditions on the gluon propagator on the 
lattice (or on the torus) we will consider twisted boundary conditions. 

As usual we denote link variables by U^{x), the unit vector in -direction by /t and the ex- 
tension of the lattice in jLt -direction by L^. Gauge fixing can be implemented by maximizing some 
functional F[g,U] with respect to gauge transformations U^{x) U^{x) = g^^ {x)Un{x)g{x + fi), 
e.g. for Landau gauge we have F[g,U] = Y,x,fi U^{x). Non-trivial boundary conditions on the lattice 
(as well as on the torus) can be introduced by demanding that the link variable is periodic up to a 
gauge transformation, i.e. 

u^, {x + Lv v) = t/J" {x) = a- ' {x)U^ {x)a^ (x+fi), (2. 1) 

where D.y{x) is called transition function. This property of the link variable ensures that gauge 
invariant observables ai^e actually periodic. The transition functions have to fulfill the cocycle 
condition 

^y{x)a^{x + LyV)=a^{x)ay{x + L^fi)Z^y, (2.2) 

where Z^y are elements of the center of the gauge group. Twisted boundary conditions imply 
that some of the group elements Z^y are non-trivial. Under a general gauge transformation g the 
transition functions transform as 

D.f,{x)=g-\x)a^,ix)g{x + L^fL) (2.3) 

leaving the twists Z^v unchanged. For twisted boundary conditions gauge fixing is non-trivial: If 
we implement non-trivial boundary conditions via some transition functions Q.^ {x) we should keep 
these transition functions fixed during gauge fixing, i.e. we have to find the maximum of the gauge 
fixing functional F[g,U] with respect to such gauge transformations g which do not change the 
transition functions. 
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But there is a second problem concerning the choice of transition functions for a given gauge 
functional F [g, ?7] . Let us consider the Landau gauge in periodic boundary conditions. Maximizing 
F[g,U] in this case means that we try to make the deviation of the gauge fixed configuration Ufi (x) 
from the configuration U^x {x) = 1 (which corresponds to the absolute maximum of F and which 
has zero field strength) as small as possible. On the other hand one could (awkwardly) choose 
non-trivial transition functions' and keep them fixed during Landau gauge fixing. Then in general 



{x) = 1 will not fulfill these new periodicity properties ( |2.1[ ) and some other configuration (x) 
will be the maximum of the functional F[l,?7]. However, in general this configuration U^{x) 
will have non-zero and non-constant field strength. Landau gauge fixing with these (artificially) 
introduced transition functions will make the deviation of the gauge fixed configuration U^{x) 
from U^{x) as small as possible. This will obviously destroy translational invariance, which is 
physically not sensible. The lesson is that one should be careful in choosing the transition functions. 
If possible one should choose the transition functions (also in the twisted case) such that (x) = 1 



fulfills eq. (gj). 

From now on we will specialize to the gauge group SU (2). We parametrize the links as usual 
(cT^: - Pauli matrices) 

U^{x) = ul{x)l + iu'l{x)Ok, w/iGR, 7 = 0,1,2,3 (2.4) 



and choose constant transition functions, so-called twist eaters 01O|]: 

D.o = l, Q.k = iOk, k = 1,2,3 (2.5) 

corresponding to the twists 

Zok = l, k= 1,2,3, Zi2 = Z23 =Z3i = -1. (2.6) 



From eq. (2.1) one easily obtains 



ul {x + Luk) = ul (x) , (x + L,) = I '^^J \ (2.7) 

i.e. some of the colour components are anti-periodic in spatial directions. The vacuum configura- 
tion {x) = 1 fulfills these boundary conditions, which are compatible with Landau gauge. 

Now let us have a look at the zero momentum part of the gauge potential. The gauge potential 
A^{x) (corresponding to the hnks U^{x)) and and its Fourier transformation are given by 



Af,ix) = l/2{Uf,ix)-U+{x)) (2.8) 
1 , 



Mp) = ^LM^y"'' (2.9) 



With periodic boundary conditions the zero-momentum part ■.=A^{p = 0) shows large fluctu- 
ations during the simulation as can be seen in fig. [l[ These fluctuations and their amplitudes 
are clearly reduced with twisted boundary conditions. At this point let us mention two further 



' This can be achieved by a non-periodic gauge transformation g. Then, according to eq. ( 2.3 ), one has the transition 
functions = g^^ {x)g{x + Lfifl) ^ 1. 
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Figure 1: Simulation history of the zero momentum part in periodic (black) and twisted (red) boundary 
conditions on a 6^ lattice with j3 = 10 (left) and j3 =2.15 (right). 
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Figure 2: Plot of gluon propagator D{p'^) with twisted and periodic boundary conditions for different 
lattice sizes (left) and comparison with results from Dyson-Schwinger equation on the Torus 0. 



advantages of the constant transition functions: They are numerically easy to implement and the 
periodicity in time allows finite temperature calculations. 

The effect of the twisted boundary conditions should become smaller on larger lattices, while 
for small lattices we expect a strong dependence on the lattice size. 

3. The gluon propagator in twisted boundary conditions 

We have measured the gluon propagator for different lattice sizes at /3 = 2.15, see fig.§ As can 
be seen: for larger spatial extension of the lattice the "twisted" gluon propagator becomes larger and 
the "twisted" gluon propagator is always smaller than the "periodic" gluon propagator, i.e. one can 
consider the gluon propagator in the twisted ensemble as a lower bound for the gluon propagator 
in the periodic ensemble. As expected for larger spatial extensions the difference between the 
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propagators with periodic and twisted boundary conditions decreases, i.e. the finite size effects due 
to the twists become smaller. 

All these observations suggest that the gluon propagator D{p^) on the lattice has a non-zero 
limit for p 0. Furthermore, our results confirm results obtained by solving Dyson-Schwinger 
equations on the torus [^. 

Unfortunately, our investigations do not clarify the difference between compact and non- 
compact space-times. Our findings seem to indicate that this difference in the infrared behaviour 
of the gluon propagator remains even in the limit of an infinitely large lattice. 
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